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ABSTRACT 


An  investigation  was  performed  to  study  the  impact  damage  resistance  of  lam¬ 
inated  composites  caused  by  a  low-velocity  foreign  object.  The  effects  of  material 
interfacial  strength  and  impactor’s  mass  and  nose  size  on  the  damage  resistance  were 
the  major  concern  of  the  study.  An  analytical  model  was  developed  for  predicting 
the  dynamic  response  of  composites  during  impact.  The  model  consists  of  damage 
accumulation  prediction  and  material  degradation  modeling.  Damage  accumulation 
criteria  were  adopted  for  predicting  failure  and  the  mode  of  failure  due  to  impact. 
Appropriate  material  constitutive  relations  were  also  proposed  to  relate  strains  to 
stresses  for  the  material  that  was  inflicted  with  impact  damage.  The  Hertzian  con¬ 
tact  law  was  also  modified  by  taking  into  account  material  degradation  during  impact. 

The  proposed  model  and  the  modified  Hertzian  contact  law  were  implemented  in 
a  transient  dynamic  finite  element  analysis,  designated  as  “3DIMPACT.”  The  code 
can  predict  the  force  response  and  damage  accumulation  of  a  composite  during  im¬ 
pact.  Composites  with  and  without  toughening  interleaves  were  tested  experimentally 
to  verify  the  model  and  the  computer  code.  The  effect  of  the  impactor’s  size  and  nose 
radius  was  also  considered  in  verifying  the  code.  Overall,  the  prediction  agreed  very 
well  with  the  test  data.  The  code  could  be  very  useful  for  evaluating  the  impact 
resistance  of  a  composite  and  also  for  screening  materials. 
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I.  INTRODUCTION 


One  of  the  major  concerns  in  the  design  of  composite  structures  is  impact  dam¬ 
age.  Organic  matrix  fiber-reinforced  laminated  composites  are  very  susceptible  to 
transverse  impaet,  esp>ecially  at  low  velocities.  Low- velocity  impact  can  cause  dam¬ 
age,  including  matrix  cracks,  delaminations,  and  fiber  breakage,  which  is  embedded 
inside  the  composites.  Such  damage  is  very  difficult  to  detect  with  the  naked  eye  and 
can  cause  significant  reductions  in  the  strength  and  stiffness  of  the  materials. 

Hence,  the  knowledge  of  damage  resistance  of  composites  is  critically  important 
for  their  appUcation.  It  has  been  observed  [1-5]  that  the  impact  damage  strongly 
depended  upon  many  factors  such  as  material  properties  of  composites,  laminate 
ply  orientation  and  thickness,  as  well  as  the  mass  and  the  shape  of  the  projectile. 
Interleaved  composites  or  composites  with  toughened  interfaces  have  been  shown 
experimentally  to  produce  higher  damage  resistance  to  impact  [6,7],  Therefore,  it 
is  possible  to  improve  the  impact  resistance  of  a  laminated  composite  by  properly 
designing  the  material  and  the  structures. 

Numerous  investigations  have  been  conducted  on  this  subject,  including  experi¬ 
mental  and  analytical  work.  Because  of  the  overwhelming  amount  of  literature  related 
to  the  subject,  only  some  of  the  studies  are  selectively  cited  here  [1-43].  Because  the 
damage  was  mostly  embedded  inside  the  materieils,  tested  specimens  were  usually 
C-scanned  or  X-rayed,  and  then  sliced  in  pieces  to  inspect  the  damage. 

Several  analytical  models  [23-43]  have  been  developed  to  study  the  transient 
dynamic  response  of  composites  due  to  impact.  However,  the  major  attention  of 
most  of  the  analyticed  studies  was  focused  either  on  the  contact  force  of  an  impactor 
or  on  the  response  of  the  plates  without  consideration  of  impact  damage.  A  recent 
study  by  the  P.I.  and  his  associates  [42,43]  has  resulted  in  a  semi-empirical  model 
as  well  as  a  computer  code,  “3D IMPACT,”  based  on  the  observed  impact  damage 
mechanisms  for  estimating  the  extent  of  impact  damage  in  laminated  composites.  No 
material  degradation  due  to  impact  was  considered  in  the  model.  By  appropriately 
selecting  an  empirical  parameter,  the  code  was  shown  to  produce  a  good  correlation 
with  the  test  data. 
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Since  impact  can  result  in  significant  damage,  which  can  degrade  material  propK 
erties  and  alter  the  impact  response  of  the  composites,  it  is  believed  that  in  order  to 
accurately  evaluate  the  impact  resistance  of  composites,  the  effect  of  damage  must 
be  considered  in  the  analysis.  Therefore,  the  objective  of  this  study  was  to  enhance 
the  ability  of  the  previous  model  by  taking  into  account  the  effect  of  damage  accu¬ 
mulation  on  the  material  and  on  the  response  of  the  structure  during  impact.  The 
model  would  be  then  implemented  in  the  computer  code  “3DIMPACT”  so  as  to  be 
able  to  assess  the  impact  resistance  of  composites  as  a  function  of  material  properties 
and  design  parameters. 

This  report  summarizes  the  results  of  the  study.  Both  the  model  and  the  exper¬ 
iments  will  be  described  briefly  and  the  comparisons  between  the  predictions  and  the 
data  will  also  be  selectively  presented.  A  detailed  description  of  the  study  and  the 
results  was  documented  in  Reference  44  as  a  Ph.D.  dissertation. 

II.  PROBLEM  STATEMENT 

Consider  a  fiber-reinforced  laminated  composite  panel  with  and  without  tough¬ 
ening  interleaves,  and  subjected  to  transverse  impact  by  a  low-velocity  spherical-nose 
impactor  as  shown  in  Figure  1.  The  ply  orientation  of  the  laminate  can  be  arbitrary 
but  must  be  symmetric  with  respect  to  its  middle  plane.  The  impactor’s  radius  and 
mass  may  vary.  For  a  given  mass  and  nose  radius  of  the  impactor,  the  following  are 
to  be  determined: 

(1)  The  velocity  of  the  impactor  required  to  initiate  the  impact  damage. 

(2)  The  force-time  history  of  the  impact. 

(3)  The  extent  of  damage  and  failure  modes  inside  the  laminate. 

(4)  The  residual  stiffness  of  the  composite  after  impact. 


III.  THE  MODEL 

The  model  consists  of  damage  accumulation  prediction  and  material  degradation 
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modeling.  Failure  criteria  were  adopted  for  predicting  the  type  and  the  extent  of  the 
damage  in  composites  due  to  impact,  and  appropriate  material  constitutive  relations 
were  proposed  to  estimate  the  residual  properties  of  damaged  composites  due  to 
impact. 


3.1  DAMAGE  ACCUMULATION  PREDICTION 

Low- velocity  impact  damage  may  consist  of  matrix  cracks,  delaminations,  and 
fiber  breakage.  It  has  been  shown  [2,3]  that  matrix  cracking  is  the  initial  damage 
mode.  Matrix  cracks  can  induce  delaminations  immediately  along  the  bottom  or 
upper  interface  of  the  cracked  layer,  depending  on  the  position  of  the  layer  in  the 
laminate.  Fiber  breakage  may  also  occur  inside  and  on  the  surface  of  the  laminate. 

Therefore,  in  order  to  predict  the  impact  damage,  a  model  must  be  capable  of 
predicting  these  three  failme  modes  in  a  correct  sequence.  Accordingly,  three  failure 
criteria  were  adopted  for  predicting  matrix  cracking,  impact-induced  delamination, 
and  fiber  breakage. 


Matrix  Cracking  Criterion 

In  order  to  predict  the  occurrence  of  the  critical  matrix  cracks,  the  matrix  failure 
criterion  proposed  previously  by  the  authors  [42,45]  will  be  adopted  here;  the  criterion 
can  be  expressed  as 


f  CM  >  1 

\  CM  <  1 
fY^Yt 
\y  =  Yc 


Failure 
No  failure 

if^j/y  >  0 
if  ffyy  ^  0 


(1) 


where  the  subscript  of  x  and  y  are  the  local  coordinates  of  the  nth  layer  parallel  and 
normal  to  the  fiber  directions,  respectively,  and  z  is  the  out-of-plane  direction,  is 
the  effective  principal  stress  in  the  ply  and  is  defined  as  follows: 
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(2) 


Yt  and  Yc  are  the  intraply  transverse  tensile  and  compressive  strengths,  respec¬ 
tively  [45],  and  S  is  the  intraply  shear  strength.  Both  Vi  and  S  strongly  depend  upon 
the  thickness  of  the  ply  group  and  the  ply  orientations  of  the  immediate  neighboring 
plies.  Vi  and  S  are  the  stresses  corresponding  to  the  first  matrix  cracking  and  can  be 
determined  based  on  the  theory  of  elasticity  and  the  principle  of  the  conservation  of 
energy.  Details  of  the  derivation  are  given  in  Ref.  [44]. 

Whenever  the  calculated  effective  stresses  in  any  one  of  the  plies  in  the  leuninate 
first  satisfy  the  criterion  (cm  =  1)  during  impact,  initial  impact  damage  is  predicted. 
It  was  assumed  then  that  the  matrix  cradc  would  propagate  throughout  the  thickness 
of  the  ply  group  which  contains  the  cracked  ply.  The  time  t  corresponding  to  the 
initial  damage  is  designated  as  tjif  •  A  delamination  could  immediately  follow  from 
the  location  of  the  matrix  crack  along  the  interfaces  of  the  ply  group.  As  the  time 
increases  (t  >  tu)  during  impact,  additional  matrix  cracking  could  be  produced  in 
the  other  layers.  Hence,  the  criterion  should  continuously  be  applied  to  the  other 
layers  for  determining  any  additional  matrix  failure.  If  no  additional  matrix  cracking 
is  found  at  any  other  layer  during  impact,  then  the  impactor’s  velocity  associated 
with  the  first  matrix  cracking  is  referred  to  as  the  impact  velocity  threshold,  which 
is  the  velocity  required  to  just  cause  the  initial  impact  damage  of  the  laminate. 


Impact-Induced  Delamination  Criterion 

It  was  believed  that  delamination  growth  due  to  low-velocity  impact  would  occvn 
only  when  the  following  two  sequential  conditions  were  met: 

(1)  one  of  the  ply  groups  immediately  above  or  below  the  concerned  interface  has 
failed  due  to  matrix  cracking  and 

(2)  the  combined  stresses  governing  the  delamination  growth  mechanisms  through 
the  thicknesses  of  the  upper  and  lower  ply  groups  of  the  interface  reach  a  critical 
value. 
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Based  on  the  above  second  statement,  the  following  impact-induced  delamination 
criterion  originally  proposed  by  Choi  and  Chang  [42]  was  adopted  and  modified  as 
follows: 


where  ”<Tyz  and  are  the  interfacial  shear  stresses  on  the  surfaces  above  and 

below  the  ply  imder  consideration,  respectively.  and  Sxz  are  the  interlaminar 
shear  strengths  in  the  direction  perpendicular  or  parallel  to  the  fiber  direction  of  the 
ply  under  consideration,  respectively.  Similar  to  Yt  and  S,  the  values  of  Syz  and  Sxz 
may  also  depend  upon  the  ply  orientation  of  the  neighboring  plies.  Determination  of 
these  distributions  as  a  function  of  the  ply  orientation  will  be  discussed  in  Section 
IV. 

Note  that  it  was  found  that  the  empirical  parameter  Da  originally  proposed  for 
Eq.  (3)  was  not  necessary  once  the  interlaminar  strength  distributions  were  deter¬ 
mined  from  experiments. 


Fiber  Breakage  Criterion 

Fiber  breakage  in  the  laminate  was  predicted  by  the  following  criterion: 


(4) 


where  Xt  is  the  longitudinal  tensile  strength  of  the  unidirectional  composite. 

Based  on  the  previous  studies  [2,3,42],  the  matrix  cracking  criterion  should  be 
apphed  first  to  predict  matrix  failime.  Once  a  critical  matrix  crack  is  predicted  in  a 
layer,  the  delamination  criteron  and  the  fiber  breakage  criterion  will  then  be  applied 
to  estimate  the  extent  of  the  delamination  and  the  fiber  breakage  areas,  respectively. 
The  local  damage  may  cause  material  degradation  in  the  damaged  ply  and  affects 
the  subsequent  dynamic  response  of  the  composite.  Therefore,  appropriate  material 
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constitutive  relations  axe  needed  to  assess  the  residual  properties  of  the  composite 
containing  impact  damage.  The  relationships  between  the  mechanical  properties  of 
the  composite  and  the  type  and  extent  of  the  damage  will  be  established  in  the  next 
section. 


3.2  MATERIAL  DEGRADATION  MODELING 


Matrix  Cracking 

Matrix  crack  density  has  been  utilized  widely  as  a  means  for  characterizing  the 
residual  properties  of  composites  containing  matrix  crack  damage.  The  constitutive 
relations  of  the  m-th  ply  in  a  laminate  as  a  function  of  crack  density  (f>  have  been 
developed  recently  by  Shahid  and  the  P.I.  [45]  and  can  be  expressed  on  the  material 
axes  as  follows: 


{Q{4>)\m  — 


/  Qxxi<f>)  Qxy{<f>)  0  \ 

QyxifP)  Qyy{4>)  0 

\  0  0  Qss{(i>)/ 


m 


(5) 


where  Qiji<f>)  =  x,y,z)  are  the  components  of  the  effective  residual  stiffness  of  a 
ply  at  a  given  crack  density  0. 

Eq.  (5)  may  vary  from  layer  to  layer,  depending  upon  the  ply  orientation  of  the 
laminate.  It  was  assumed  in  Eq.  (5)  that  the  out-of-plane  engineering  constants,  Ezz, 
Gxz,  Gyzz,  i^xz,  and  ^yz,  were  imaffected  by  matrix  cracking. 

The  crack  density  in  Eq.  (5)  would  increase  if  the  stresses  continue  to  increase. 
Because  matrix  cracks  accumulate  rapidly  during  impact,  it  was  assumed  that  once 
the  matrix  cracking  criterion  is  predicted,  cracks  would  inunediately  reach  the  satura¬ 
tion  density,  0o,  within  the  area  where  the  stresses  satisfy  the  criterion.  Accordingly, 
the  effective  engineering  properties  of  each  degraded  ply  can  be  expressed  in  terms  of 
Qij{(f>o)  of  the  effective  ply  stiffness  as  follows: 
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riD  _  Qxxi<f>o) 

■^xx 

(6.1) 

m 

Qyyi^o) 

^yy 

m 

(6.2) 

D  Qxy{4>o) 

"  “  QvuiM 

(6.2) 

D  Qxy{<f>o) 

(6.3) 

(6.4) 

where  m  =  [1  —  and  the  superscript  D  denotes  degraded  properties. 


Impact-Induced  Delamination 

If  delamination  failure  is  predicted  after  matrix  cracking  has  occurred,  cracks 
would  be  generated  along  the  interfaces  of  the  cracked  ply,  resulting  primarily  in 
reduction  of  out-of-plane  stiffness  of  the  laminate.  Since  the  delamination  is  induced 
by  intraply  matrix  cracks,  the  effect  of  the  delamination  on  the  laminate  stiffness 
can  be  talcen  into  account  by  appropriately  reducing  the  out-of-plane  stiffness  of  the 
cracked  ply. 

Once  a  ply  has  suffered  both  interply  and  intraply  cracks,  the  actual  effective 
stiffness  of  the  cracked  ply  in  the  composite  is  very  difficult  to  be  determined  anedyt- 
ically  or  experimentally.  Since  the  cracks  have  saturated  in  the  ply,  it  was  assumed 
in  the  proposed  model  that  such  a  cracked  ply  in  the  damaged  area  can  be  treated  as 
a  transversely  isotropic  medivun.  As  a  result,  the  out-of-plane  engineering  properties 
of  the  cracked  ply  could  have  the  following  simple  expressions: 


E: 


D 


Qxx{4^o) 

m 

Qyyi^o) 

m 

Qxy{4>o) 

Qyy{4>o) 


(7.1) 

(7.2) 

(7.3) 


7 


G. 

E: 


yx 

D 

xy 

D 

zz 


V, 


xz 
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xz 


U. 


D 


zx 


yx 


Qxyi<l>o) 

Qxxi^o) 

(7.4) 

—  Qsa{4^o) 

(7.5) 

_  jpD 

-^yy 

(7.6) 

=  Gg 

(7.7) 

=  ^xy 

(7.8) 

(7.9) 

=  i'zy 

(7.10) 

where  no  degradation  was  assumed  for  Uyz  for  preserving  the  assumption  of  transverse 
isotropy. 


Fiber  Breakage 

If  the  fiber  breakage  criterion  is  predicted,  the  amount  of  the  material  degra¬ 
dation  depends  upon  the  size  of  the  area  where  the  stresses  satisfy  the  fiber  failure 
criterion.  Based  on  the  concept  of  continuum  mechanics,  a  degradation  parameter  df 
will  be  multipled  with  all  the  engineering  properties  [45],  where  df  can  be  expressed 
as 


ds  =  (8) 

where  Af  is  the  area  where  the  stresses  satisfy  the  fiber  breakage  criterion,  and  6  is 
the  fiber  interaction  length  which  can  be  estimated  from  the  micromechanics  [45].  /3 
is  the  rate  of  the  material  degradation  due  to  fiber  breakage.  Any  value  greater  than 
8  was  appropriate  for  graphite/epoxy  composites. 

It  is  worth  noting  that  the  sequence  in  applying  the  failure  criteria  and  material 
constitutive  relations  is  very  important  for  accurately  predicting  impact  damage. 
Matrix  cracks  should  be  evaluated  first.  If  matrix  cracking  is  predicted,  the  impact- 
induced  delamination  and  fiber  breakage  criteria  should  then  be  applied  to  estimate 
the  extent  of  the  delamination  and  fiber  breakage  areas.  The  procedure  has  to  be 
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rep>eated  at  the  other  layers  during  impact.  The  final  size  of  each  delamination  is 
determined  by  the  area  within  which  the  stress  components  satisfy  the  delamination 
failure  criterion  during  the  entire  duration  of  impact. 


IV.  THE  HERTZIAN  CONTACT  LAW 


The  Hertzian  contact  law  has  been  widely  utilized  to  model  the  force  response 
of  an  impactor  during  impact.  The  force-indentation  relationship  for  low-velocity 
impact  (no  penetration)  is  normally  generated  by  a  static  indentation  test.  The 
result  of  the  test  data  was  then  fitted  to  the  contact  law  in  order  to  select  a  power 
constant  for  a  particular  target  material.  However,  experiments  have  shown  that  the 
force-indentation  relationship  upon  loading  is  different  from  unloading.  In  general, 
upon  loading,  the  Hertzian  contact  law  can  be  expressed  [42,46]  as 


F  =  Kof  (9) 

Thus,  the  contact  force  F  can  be  related  to  the  indentation  depth  a  (the  distance 
between  the  center  of  the  projectile’s  nose  and  the  mid-surface  of  the  plate),  k  is  the 
contact  stiffness,  and  7  is  the  power  constant  which  has  to  be  determined  from  the 
static  identation  test. 

The  contact  stiffness  k  is  given  as 


-  4  /-  1 

''  3^  [ks  +  kt] 


(10) 


where  r  is  the  radius  of  the  impactor.  kt  can  be  related  to  the  properties  of  the  target, 
while  ks  is  related  to  projectile  properties  and  is  given  as 


ka  — 


(1  - 

E, 


(11) 


where  Vg  and  Eg  are  the  Poisson’s  ratio  and  the  Young’s  modulus  of  the  impactor. 


9 


respectively. 

For  laminated  composites,  kt  can  be  expressed  directly  in  terms  of  laminate 
properties  [47, 48].  However,  once  damage  occurs  in  the  laminate,  both  kt  and  7  may 
not  remain  as  constants,  especially  for  kt,  which  may  be  a  function  of  the  type  and 
the  extent  of  the  damage  in  the  target. 

Accordingly,  in  this  study,  the  value  of  kt  was  recalculated  at  every  time  step 
based  on  the  degraded  material  properties.  The  detailed  derivation  and  implementa¬ 
tion  of  the  kt  in  the  analysis  is  described  in  Reference  44. 

Upon  unloading,  the  predicted  contact  force  based  on  the  contact  law  in  Elq. 
(9)  deviates  from  the  measured  contact  force  for  a  given  indentation.  It  has  been 
shown  [49]  that  the  contact  law  in  Eq.  (9)  produces  a  longer  impact  diuation  upon 
imloading.  Based  on  the  test  data,  an  empirical  expression  was  suggested  by  Sun 
et  al  [32,33]  for  the  force-identation  relation.  However,  such  an  empirical  expression 
could  not  be  generalized  and  is  difficult  to  apply  once  damage  starts  to  accumulate 
in  the  target. 

Therefore,  in  the  study,  the  contact  law  in  Eq.  (9)  was  modified  to  take  into 
accoimt  local  damage,  but  was  still  used  for  both  loading  and  unloading  cases.  Accord¬ 
ingly,  it  is  expected  that  using  the  modified  contact  law,  the  dynamic  force  response 
upon  unloading  may  not  be  accurate.  This  may  affect  the  accuracy  of  the  stresses 
and  deformations  of  the  target  upon  unloading.  However,  it  will  be  demonstrated 
that  most  impact  damage  is  induced  upon  loading.  Since  the  p>ost-impact  dynamic 
response  is  out  of  the  scope  of  the  interest,  this  approach  was  adequate  for  the  present 
study. 


V.  THE  FINITE  ELEMENT  ANALYSIS 


The  “3D IMPACT”  computer  code,  previously  developed  by  Choi  and  Chang  [42], 
was  modified  by  implementing  the  proposed  model.  The  finite  element  analysis  was 
needed  for  calculating  the  stresses  and  strains  inside  the  composites  during  impact. 
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The  information  regarding  finite  element  procedures  was  given  extensively  in  [35,42]. 
Here,  only  a  brief  description  of  the  analytical  approach  will  be  given  as  follows: 

The  analysis  was  based  on  a  three-dimensional  linear  elasticity  theory.  The 
materials  in  each  layer  were  considered  homogeneous  and  orthotropic.  Accordingly, 
the  equihbriiun  equations  at  instant  time  t  in  a  variational  form  can  be  expressed  as 
[42] 


0—  Wipui^ttdv  +  /  eijEijkimdv  -  /  WiaijUjdA  (12) 

JQ  Jn  Jr 

where  aij  are  the  stresses,  eju  are  the  strains,  p  is  the  density,  Ui^tt  are  the  accelerations 

=  d^Uildt^),  Wi  are  the  arbitrary  variational  displacements,  eij  are  the  strains 
from  the  arbitrary  variational  displacements,  is  the  entire  plate  volmne,  T  is  the 
surface  of  the  plate,  Uj  is  the  outward  unit  normal  vector  on  the  plate  surface,  and 
Eijki  are  the  material  properties  of  the  laminate,  which  vary  from  layer  to  layer 
according  to  the  ply  orientation  of  the  composite. 

It  has  been  shown  that  thermal  residual  stresses  could  reduce  the  impact  resis¬ 
tance  of  composites  [2,3].  The  thermal  stresses  due  to  fabrication  were  calculated 
based  on  the  theory  of  hnear  elasticity.  The  thermal  stresses  were  combined  with  the 
mechanical  stresses  calculated  from  Eq.  (12)  to  determine  the  damage  state  in  the 
failure  criteria. 

In  order  to  solve  Eq.(l),  the  distribution  of  the  contact  force,  F{=  orijUj),  be¬ 
tween  the  impactor  and  the  impacted  laminate  must  first  be  known.  The  Hertzian 
contact  law  was  utilized  to  generate  the  contact  force.  However,  the  contact  stiffness 
was  adjusted  according  to  the  amoimt  of  damage  inside  the  target. 

An  eight-node  brick  element  incorporating  incompatible  modes  developed  previ¬ 
ously  by  Wu  and  Chang  [35]  was  used  in  the  finite  element  calculations.  The  inclusion 
of  incompatible  modes  weis  to  improve  the  accuracy  in  calculating  the  bending  stiff¬ 
nesses  and  the  interlaminar  shear  stresses.  A  direct  Gauss  quEidrature  integration 
scheme  was  adopted  [39]  through  the  element  thickness  to  account  for  the  change 
in  material  properties  from  layer  to  layer  within  the  element.  Therefore,  plies  with 
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different  ply  orientations  could  be  grouped  in  an  element,  resulting  in  a  significant 
reduction  in  computational  time  and  memory  space  for  the  3-dimensional  analysis. 

The  contact  area  may  increase  upon  loading  and  decrease  upon  unloading  during 
the  impact.  The  contact  force  was  assumed  to  be  distributed  in  a  normal  distribution 
over  the  area  [44]. 


IV.  EXPERIMENTS 

4.1  INTERLAMINAR  SHEAR  STRENGTHS 

Experiments  were  conducted  to  measure  the  interlaminar  shear  strengths,  Sxz 
and  Syz,  as  a  function  of  the  ply  orientation  mismatch.  Information  about  these  inter¬ 
laminar  shear  strengths  is  required  from  the  impact-induced  delamination  criterion. 
The  shear  strengths  Sxz  and  Syz  of  a  ply  in  a  laminate  are  defined  as  the  interfacial 
shear  strengths  along  the  fiber  direction  or  normal  to  the  fiber  direction  of  the  ply 
under  consideration,  respectively.  Accordingly,  depending  upon  the  ply  orientation, 
the  interlaminar  shear  strengths  may  vary  with  the  ply  orientation  mismatch. 

The  short  beam  shear  test  was  selected  for  measming  the  interlaminar  shear 
strength.  In  order  to  evaluate  the  effect  of  ply  angle  mismatch  on  the  shear  strengths, 
the  ply  orientations  of  the  specimens  were  carefully  selected  for  the  tests. 

The  interlaminar  shear  strength  was  taken  as  the  stress  at  the  interface  corre¬ 
sponding  to  the  collapse  load  of  the  specimens.  In  order  to  calculate  shear  stresses  on 
the  interfaces,  the  classical  lamination  theory  was  adopted  to  calculate  the  inplane 
stresses  based  on  the  three-point  bending  solutions.  The  out-of-plane  shear  stresses 
were  then  determined  by  employing  the  three-dimensional  equilibrium  equations  with 
appropriate  boundary  conditions. 

A  summary  of  the  interlaminar  shear  strength  distributions  between  two  neigh¬ 
boring  plies  as  a  function  of  the  mismatch  of  the  ply  angles  is  shown  in  Figures  2  and 
3.  Apparently,  the  difference  in  ply  angles  affects  significantly  the  interlaminar  shear 
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strength.  Additionally,  because  the  distributions  of  Syz  and  Sxz  are  considerably  dif¬ 
ferent,  it  is  expected  that  the  growth  of  interfacial  delamination  can  be  different  in 
the  direction  parallel  and  normal  to  the  fiber  direction  of  the  ply  under  consideration. 


4.2  COMPOSITES  WITH  AND  WITHOUT  INTERLEAVES 

In  order  to  verify  the  proposed  model,  composites  with  and  without  toughening 
interleaves  were  fabricated.  T300/976  composites  and  FM300  thermoset  interleaves 
were  selected  for  the  study.  The  thickness  of  a  single  interleaf  is  about  one  fourth 
of  the  composite  prepreg  thickness.  Numerous  ply  orientations  were  chosen.  The 
dimensions  of  each  specimen  were  10  cm  long  and  7.6  cm  wide. 

All  the  specimens  were  cured  without  thermal-induced  pre-matrix  cracks.  AU 
the  specimens  were  cut  by  a  diamond-coated  saw  and  X-rayed  after  the  cutting  to 
inspect  any  internal  damage  due  to  manufacturing  and  cutting.  No  apparent  damage 
was  found  from  the  X-radiographs  after  cutting  all  the  specimens. 

Figure  4  shows  schematically  the  impact  test  facility  used  for  the  study.  A 
spherical-nosed  impactor  was  selected  for  the  study.  The  radius  of  the  spherical  nose 
head  made  of  steel  was  0.635  cm.  The  specimens  were  firmly  clamped  along  two 
parallel  edges  as  shown  in  Figure  4. 

The  results  of  the  X-radiograph  examinations  of  all  the  tested  specimens  showed 
consistently  that  the  laminates  with  interleaves  sustained  higher  energy  impact  at 
the  same  damage  size.  The  selected  results  of  the  comparison  are  shown  in  Figures 
5-7.  It  is  noted  that  the  double  slashes  “//”  appearing  in  the  ply  orientations  of 
these  figures  represented  an  interleaf  at  the  location.  The  damage  area  in  the  figures 
was  estimated  directly  from  the  X-radiographs.  Apparently,  both  impact  velocity 
threshold  and  impact  damage  size  were  sensitive  to  the  presence  of  the  interleaves. 
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V.  COMPARISON  AND  VERIFICATION 


5.1  EFFECT  OF  DAMAGE  ON  FORCE  RESPONSE 

Test  data  available  in  the  literature  were  selected  and  utilized  to  verify  the  model 
and  the  computer  code.  First,  the  predictions  were  compared  with  the  data  generated 
from  heavy  mass /low  velocity  impact.  Secondly,  the  effect  of  the  impactor’s  mass  on 
the  impact  response  was  evaluated  for  a  given  impact  energy. 


Heavy  Mass-Low  Velocity  Impact 

Figures  8-10  compared  the  force-time  history  of  a  heavy  impactor  on  IM7/954- 
2  composites  between  the  predictions  and  the  measurements.  The  test  data  were 
obtained  from  Lockheed  [51].  Numerical  simulations  of  the  test  condition  were  per¬ 
formed  with  and  without  the  proposed  impact  damage  model.  The  results  of  both 
calculations  were  presented  in  the  figures  and  compared  with  the  test  data. 

Upon  the  loading  path,  the  predictions  based  on  the  present  model  agreed  with 
the  data  very  well.  The  first  load  drop  from  the  data  corresponded  quite  well  with 
the  predicted  load  at  which  the  initial  damage  was  predicted  due  to  matrix  cracking 
and  delamination.  The  calculations  without  considering  the  damage  could  not  detect 
the  load  drop  and  overestimated  the  maximum  load  that  the  composite  could  sustain. 

Upon  unloading,  both  predictions  with  and  without  the  inclusion  of  the  damage 
model  deviated  from  the  data.  As  expected,  the  deviation  can  be  attributed  primarily 
to  the  Hertzian  contact  law,  which  was  developed  originally  not  for  simulating  the 
unloading  response.  However,  by  examining  the  predicted  size  of  the  accumulated 
damage,  it  was  foimd  that  most  of  the  damage  in  composites  was  generated  along 
the  loading  path.  Unloading  could  affect  the  post-impact  dynamic  response  of  the 
target,  which  is  not  of  particular  interest  to  the  study. 

Figure  1 1  shows  the  contact  stiffness  distribution  of  the  impactor  as  a  function 
of  time  during  impact  at  three  different  velocities.  Clearly,  the  contact  stiffness 


14 


degraded  once  damage  occurred  in  the  laminate.  At  higher  velocity  at  which  damage 
accumulated  faster,  the  rate  of  degradation  was  high.  Upon  miloading,  the  contact 
stiffness  remained  almost  constant  for  all  the  cases  studied. 


Mass  Effect 

The  effect  of  mass  on  the  impact  response  of  T800/3900-2  composites  at  a  given 
impact  energy  was  studied  by  Delfosse  et  al.  [49].  Figures  12-15  showed  the  predicted 
force-time  response  as  compared  to  the  measurements  taken  from  [49].  At  the  same 
impact  energy,  the  impactor’s  mass  affected  significantly  the  force  response.  At  a 
heavier  mass,  a  single  hump  was  foimd,  but  multiple  humps  were  found  for  a  lower 
mass.  Again,  the  predictions  matched  with  the  data  very  well  on  the  loading  path. 
Upon  unloading,  the  model  predicted  a  longer  duration  of  impact  than  was  recorded 
in  the  experiments. 


5.2  TOUGHENED  V.S.  UNTOUGHENED  COMPOSITES 

Numerical  simulations  were  also  generated  to  compare  with  the  impact  test  data 
that  were  generated  on  composites  with  eind  without  interleaves  dmring  the  inves¬ 
tigation.  Figures  5-7  also  show  the  comparison  between  the  predictions  and  the 
measurements  of  the  damage  size  as  a  function  of  impact  energy.  Clearly,  the  pre¬ 
dictions  agreed  quite  well  with  the  measurements  for  the  laminates  with  and  without 
interleaves. 


VI.  COMPUTER  CODE:  3DIMPACT 


The  computer  code,  “3DIMPACT,”  has  been  considerably  modified  based  on  the 
proposed  model.  The  code  has  a  user-firiendly  interactive  mode  for  inputs.  The 
code  was  developed  for  analyzing  laminated  composite  plates  subjected  to  transverse 
projectile  impact.  For  a  given  mass  and  raduis  of  the  impactor  and  the  material 
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properties  of  a  composite,  the  code  is  capable  of  producing  the  following  results: 

1  )  three-dimensional  stress  and  strain  distributions  during  the  impact. 

2  )  the  impact  force-time  history. 

3  )  the  mode  and  the  extent  of  the  damage  in  every  layer  of  the  laminate  due  to 
impact. 

4  )  the  residual  stiffness  of  the  composite  after  impact. 

A  copy  of  the  code  can  be  obtained  from  Pu-Kuo  Chang  at  the  address  given  on 
the  cover  page. 


VII.  CONCLUSIONS 


A  model  has  been  developed  for  predicting  impact  damage  and  residual  prop¬ 
erties  of  laminated  composites.  The  model  has  been  implemented  in  “3D IMPACT” 
code.  Experiments  were  also  conducted  on  laminates  with  and  without  toughening 
interleaves  to  verify  the  model  and  the  computer  simulations. 

The  predictions  based  on  the  model  agreed  with  the  measured  force  response 
upon  loading  very  consistently.  However,  the  predictions  overestimated  the  dmation 
of  the  impact  upon  imloading,  which  could  be  attributed  to  the  Hertzian  contact  law. 
Meanwhile,  the  model  predicted  very  well  the  impact  damage  size  for  materials  with 
and  without  interleaves. 

Based  on  the  study  the  following  remarks  can  be  made: 

(1)  The  force-time  response  of  large-mass  impact  is  more  sensitive  to  impact  damage 
than  low-mass  impact  at  the  same  impact  energy. 

(2)  Alternating  ply  orientations  produce  smaller  damage  size  than  clustered  ply 
orientations  of  the  same  thickness. 

(3)  Composites  with  toughened  interfaces  result  in  higher  impact  resistance. 
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(4)  Interfacial  shear  strengths  as  well  as  transverse  tensile  strength  are  critical  in 
governing  impact  damage  size. 

It  is  believed  that  the  Hertzian  contact  law  has  to  be  further  modified  in  order 
to  study  the  unloading  path  during  impact.  A  direct  modelling  of  the  impactor  as  an 
elastic  or  rigid  body  in  the  finite  element  simulation  would  be  desired  to  avoid  the 
uncertainty  in  the  contact  law. 
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Figure  1  Description  of  the  problem. 
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Figure  2  Interlaminax  shear  strength  Sxz  of  a  unidirectional  ply  in  a  laminate  as  a  function 
of  the  ply  angle  mismatch  with  respect  to  its  neighboring  ply. 


Material:  T300/976 
Interleaf:  FM300 


20 


Figure  3  Interlaminar  shear  strength  Syz  of  a  unidirectional  ply  in  a  laminate  as  a  function 
of  the  ply  angle  mismatch  with  respect  to  its  neighboring  ply. 
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Schematic  of  the  impact  test  set-up. 
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act  damage  area  as  a  function  of  the  impact  velocity  for  composites  with 
without  interleaves.  Comparison  between  the  predictions  and  the  test  data. 


(^inui)  vaav  aaDVPMva 


23 


Figure  6  Impact  damage  area  as  a  function  of  the  impact  velocity  for  composites  with 
and  without  interleaves.  Comparison  between  the  predictions  and  the  test  data. 
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Figure  7  Impact  damage  area  as  a  function  of  the  impact  velocity  for  composites  with 
and  without  interleaves.  Comparison  between  the  predictions  and  the  test  data. 
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Figure  9  Impact  force  distribution  as  a  function  of  time  for  IM7/954-2  graphite /epoxy 
composites.  Comparison  between  the  predictions  with  and  without  material 
degradation  and  the  test  data. 
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Figure  10  Impact  force  distribution  as  a  function  of  time  for  IM7/954-2  graphite/epoxy 
composites.  Comparison  between  the  predictions  with  and  without  materied 
degradation  and  the  test  data. 
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Figure  13  Impact  force  distribution  as  a  function  of  time  for  T800/3900-2  graphite/epoxy 
composites.  Comparison  between  the  predictions  based  on  the  present  model 
and  the  test  data  [491. 
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Figure  14  Impact  force  distribution  as  a  function  of  time  for  T800/3900-2  graphite /epoxy 
composites.  Comparison  between  the  predictions  based  on  the  present  model 
and  the  test  data  [49]. 
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Figure  15  Impact  force  distribution  as  a  function  of  time  for  T800/3900-2  graphite/epoxy 
composites.  Comparison  between  the  predictions  based  on  the  present  model 
and  the  test  data  [49]. 
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